Superexchange Interactions in Orthorhombically 
Distorted Titanates OTiOg {R= Y, Gd, Sm, and La) 

I V Solovyev 

Computational Materials Science Center, National Institute for Materials Science, 
1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan 

E-mail: solovyev.igor@iiims.go.jp 

Abstract. Starting from the multiorbital Hubbard model for the t2g-bands of i?Ti03 
{R— Y, Gd, Sm, and La), where all the parameters have been derived from the first- 
principles electronic structure calculations, we construct an effective superexchange 
(SE) spin model, by treating transfer integrals as a perturbation. We consider four 
approximations for the SE interactions: (i) the canonical crystal-field (CF) theory, 
where the form of the the occupied t2£(-orbitals is dictated by the CF splitting at 
each Ti-site and three extensions of the CF theory, namely (ii) the relativistic one, 
where occupied orbitals are confined within the lowest Kramers doublet obtained from 
the diagonalization of the crystal field and relativistic spin-orbit (SO) interactions; 

(iii) the finite-temperature extension, which consider the effect of thermal orbital 
fiuctuations near the CF configuration on interatomic interactions between the spins; 

(iv) the many-electron extension, which is based on the diagonalization of the full 
Hamiltonian constructed in the basis of two-electron states separately for each bond 
of the system. The main results are summarized as follows, (i) Thermal fiuctuations 
of the orbital degrees of freedom can substantially reduce the value of the magnetic 
transition temperature, (ii) The relativistic SO coupling is generally responsible for 
anisotropic and antisymmetric Dzyaloshinsky-Moriya interactions. All interactions 
are rigorously derived and their implications to the magnetic properties of i?Ti03 are 
discussed, (iii) The CF theory, although applicable for YTiOa and high-temperature 
structures of GdTiOa and SmTiOa, breaks down in the case of LaTiOa. In the latter, 
the CF splitting is small. Therefore, the many-electron effects in the bonds as well 
as the relativistic SO interaction start to play an important role. It is argued that 
the combination of these two effects can be responsible for the AFM character of 
interatomic correlations in LaTiOa. (iv) The SE interactions in YTiOa strongly depend 
on the details of the crystal structure. Crystal distortions in the low-temperature 
structure tend to weaken the ferromagnetic interactions. 



PACS numbers: 75.10.Dg, 75.30.Et, 71.70.Ej, 71.10.-w 



1. Introduction 

Recently, the titanium perovskites i^TiOs (where R is the rare-earth element or Y) 
have attracted considerable attention. These, formally isostructural and isoelectron 
materials, crystallize in the orthorhombic space group D^f^ and have only one electron 
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in the t2g-she\\ of Ti-atoms. Nevertheless, the magnetic properties of RTiO^ appear to be 
extremely sensitive to the magnitude and details of the lattice distortion. More distorted 
compounds with R= Yb - Gd (including YTiOs) appear to be ferromagnetic (FM). The 
Curie temperature (Tc) typically varies from about 30 K {R= Y and Gd) till 60 K {R= 
Dy, Ho, and Tm) [HE], [3]. Less distorted compounds {R= Sm - La) form the so-called 
G-type antiferromagnetic (AFM) structure, where all nearest-neighbor (NN) spins are 
coupled antiferromagnetically. The Neel temperature (Tn) varies monotonously from 
146 K (LaTiOs) till 48 K (SmTiOs). The complete phase diagram of i^TiOs can be 
found in [2]. 

The origin of the FM-AFM transition and, more generally, the microscopic physical 
picture underlying the behavior of these Mott-Hubbard insulators with partially filled 
t2g-shell is a matter of considerable debates. Different types of theories ranging from the 
degenerate orbital liquid [3] to the distortion-controlled crystal-field (CF) splitting of the 
atomic t2g-levels O IH [7] have been proposed. The main obstacle with the setting of the 
proper physical model is related to the fact that there are too many model parameters, 
which can drastically change the picture if they are chosen in a uncontrollable way. 
For example, the orbital liquid state is rapidly deteriorated by the lattice distortion 
away from the simple cubic structure. The conventional CF theories [Sj [7] crucially 
depend on the choice of the dielectric constant, which controls the magnitude of the 
t2g-level splitting. Thus, in unbiased theories for i^TiOa, it is very important to reduce 
to the minimum all arbitrariness related to the choice of the model parameters by using 
for these purposes first-principles electronic structure calculations. This direction is 
sometimes called "the realistic modeling of strongly correlated systems" [8]. 

Our previous works P, [9l [10] were devoted to construction of the parameter-free 
lattice fermion model (more specifically - the multiorbital Hubbard model) for the low- 
energy t2g-bands of distorted perovskite oxides: 

The model itself is specified in the basis of six Wannier orbitals (the so-called spin- 
orbitals), which are denoted by Greek symbols, each of which is a combination of two 
spin (s= t or |) and three orbital (m= 1, 2, or 3) variables. The site-diagonal part 
of hij = includes the CF Hamiltonian hf^ and the relativistic spin-orbit (SO) 

interaction h^^ . The off-diagonal {i^j) elements of hij stand for the transfer integrals 
iij = Both hf^ and iij are diagonal with respect to the spin indices. Uap-ys are the 

matrix elements of screened Coulomb interactions. All parameters were determined in 
the ab initio fashion, in the frameworks of the density functional theory. Other details 
can be found in the original work P3] as well as in the review article 0. 

In this paper we discuss how the lattice fermion model ([1]) for the t^g compounds 
can be further mapped onto the the spin- 1/2 model of the general form 

7t;5 = ^(o-,,4(t,), (2) 

(ij) 
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where Aij is the three-dimensional tensor, describing various interactions between 
spins at the sites i and j. It can be further decomposed into isotropic Heisenberg 
(Jij), antisymmetric Dzyaloshinsky-Moriya (dij), and symmetric anisotropic (fjj) 
interactionsjj so that the model itself can be rewritten 

ns = - .hj{cri, a-j) + ^(dij, [&i X aj]) + ^ia-i, Tijo-j), (3) 

(ij) (ij) (ij) 

where <Tj = {af, erf, a^) is the vector of Pauli matrices in the spin subspace and the 
prefactor S*^ for 5=1/2 is included in the definition of the model parameters Aij, Jij, dij, 
and fij. Thus, Jij is the scalar, d^j = [d^-,d\j,dlj) is the vector, and fij = Wt""^]] (where 
a and b stand to denote the x-, y-, or 2;-components in the orthorhombic coordinate 
frame) is the three-dimensional tensor satisfying the condition r,^-^ + rfj^ + t^j^ = 0. The 
notation (<Tj, <Tj) stands for the inner product and [ai x a-j] - for the cross product of 
the vectors ct, and a-j. 

The goal of this work is twofold. On the one hand we will discuss different levels of 
approximations for the superexchange (SE) interactions underlying derivation of the spin 
Hamiltonian ([2]). On the other hand, we investigate the accuracy of realistic modeling 
for different types of materials. Particularly, how far one can go with the description of 
not only qualitative but also quantitative aspects of interatomic magnetic interactions 
in i^TiOa. For these purposes we will consider four characteristic compounds: FM 
YTiOs and GdTiOs, and G-type AFM SmTiOs and LaTiOa. We will also present 
details of the parameters of the Hubbard Hamiltonian and interatomic magnetic 
interactions (EHS]). After brief introduction of the crystal structure of -RTiOs (Sec. [2]) 
and description of parameters of the Hubbard Hamiltonian as derived from the first- 
principles electronic structure calculations (Sec. [3]), we will turn directly to the analysis 
of the SE interactions. We will start with the one-electron approximation (referring to 
the form of the ground-state wavefunction in the atomic limit) and consider different 
types of the one-electron theories, such as the conventional CF theory at T=0 (Sec. 14. ip . 
CF theory with the relativistic SO interactions, which explains the appearance of 
anisotropic and antisymmetric Dzyaloshinsky-Moriya interactions (Sec. 14.21) . We will 
also estimate the effect of thermal fluctuations of the orbital degrees of freedom on the 
interatomic interactions between the spins (Sec. [5]). In Sec. IHwe will consider a many- 
electron extension of the SE theory, where the ground-state wavefunction is determined 
from the diagonalization of two-electron Hamiltonians constructed separately for each 
bond of the system. Particularly, we will argue that the behavior of LaTiOa is drastically 
differrent from other compounds. If in YTiOa, GdTiOa, and SmTiOs the type of the 
magnetic ground state is mainly determined by the CF splitting, which is clearly larger 
than the energy gain caused by the virtual hoppings underlying the SE processes and the 
relativistic SO interaction, in LaTiOs all three factors are at least comparable making 
the situation less straightforward. This revives the idea of our earlier work [9j and now 

^ Note that for the spin 1/2, the single- fon anisotropy term is absent. 

^ Note however that the CF sphtting calculated in l9j is incorrect, because it does not take into 
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we will explicitly show how the combination of many-electron effects in the bonds and the 
relativistic SO interaction may explain the AFM character of interatomic correlations 
in LaTiOs- Finally, brief summary of the work and main conclusions will be presented 
in Sec. |6l 

2. Crystal Structure 

The titanates i^TiOs {R= Y, Gd, Sm, and La) crystallize in the orthorhombic space 
group is Dlfi^ in Schonflies notations (No. 62 in International Tables). An example of 
the crystal structure with the notations of four Ti-atoms composing the primitive cell 
is shown in Fig. [T]for the case of YTiOs |12j . 




Figure 1. Crystal structure of YTiOa. The Y atoms are indicated by the big blue 
spheres, the Ti are indicated by the medium red spheres, and the oxygen atoms 
are indicated by the small green spheres. The symbols a, b, and c stand for the 
orthorhombic translations. Different Ti-atoms are shown by numbers. 

All calculations have been performed in the experimental crystal structure. The 
low-temperature (LT) structure is currently available only for YTiOa and LaTiOa, which 
was measured at T= 2 K [2j and T= 8 K [13], respectively. Moreover, in order to make a 
direct comparison with our previous works [8[ , we also used another set of parameters 
for YTiOa, corresponding to the room temperature (RT) [12]. For SmTiOa and GdTiOa, 
the crystal structure was taken from [2], correspondingly for T= 100 K and 290 K. 

Other details of the crystal structure and their implications to the properties of 
t2g-compounds can be found in the previous publications [8i ilOj . 

account the nonsphericity of the Madelung potential [5j. Similar problem exists in [6]. The situation 
was discussed in ^10) . 
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3. Parameters of Hubbard Model 

In this section we summarize parameters of the Hubbard model ([T]), as obtained from the 
first-principles electronic structure calculations. The procedure used for the construction 
of the low-energy model ([T]) as well as for the calculation of the model parameters was 
described in many details in [HI HI]. It is true that the procedure itself relies on a 
number of approximations. Nevertheless, we would like to emphasize that apart from 
these approximations, we did not use any adjustable parameters. In this sense, the 
procedure is parameter-free. 

Basically, there are three sets of the model parameters describing 

(i) the CF sphtting of the atomic t2g-levels, 

(ii) the transfer integrals, 

(iii) the on-site Coulomb and exchange interactions. 

The behavior of parameters of the crystal field is explained in table [H which shows 
the eigenvalues and eigenvectors obtained after the diagonalization of hf^. One can 

Table 1. Eigenenergies (measured in meV from the lowest energy) and eigenvectors 
obtained after the diagonaUzation of the crystal-field Hamiltonian hf^. The 
eigenvectors are expanded over the basis of the xy, yz, z'^, zx, and x^-y'^ orbitals, 
in the orthorhombic coordinate frame. 



compound 


energies 




orbitals 











-0.13, 


0.45, 


0.38, 


-0.61, 


0.51 


YTiOg (RT) 


109 


0.46, 


0.70, 


0.14, 


0.08, 


-0.52 




123 


0.06,^ 


-0.46, 


0.32, 


-0.61, 


-0.56 







-0.19, 


0.40, 


0.33, 


-0.57, 


0.61 


YTiOs (LT) 


144 


0.10,- 


-0.33, 


0.37, 


-0.66, 


-0.57 




226 


0.43, 


0.80, 


0.15, 


0.09, 


-0.39 







-0.17, 


0.43, 


0.29, 


-0.55, 


0.64 


GdTiOa 


126 


0.17,- 


-0.14, 


0.37, 


-0.67, 


-0.60 




191 


0.38, 


0.85, 


0.04, 


0.21, 


-0.30 







-0.16, 


0.50, 


0.21, 


-0.42, 


0.71 


SmTiOa 


112 


0.13,- 


-0.09, 


0.38, 


-0.77, 


-0.49 




167 


0.38, 


0.82, 


0.02, 


0.21, 


-0.37 







-0.06,- 


-0.85, 


-0.15, 


0.34, 


-0.37 


LaTiOa 


37 


0.20, 


0.46, 


-0.26, 


0.78, 


-0.28 




61 


-0.27,- 


-0.15, 


-0.17, 


0.41, 


0.84 



clearly see that the CF splitting tends to quench the orbital degrees of freedom for all 
considered compounds, except LaTiOa: the lowest level, split off by the crystal field, 
is nondegenerate and separated from the next level by an energy gap of at least 100 
meV. Thus, in the atomic limit, the single electron occupies the lowest t2g-leve\, and 
the degeneracy of the ground state is lifted by the crystal field. The CF splitting 
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decreases in the direction YTiOa— i^GdTiOa— s^SniTiOs— s^LaTiOa. Especially, in LaTiOs 
the CF splitting is exceptionally small. In YTiOs, there is a clear dependence of the 
CF splitting on the crystal structure. Basically, there are two effects: (1) the order of 
the middle and highest t2g-levels is reversed in the RT structure (in comparison with 
the LT one); (2) the CF splitting itself is substantially smaller in the RT structure. The 
structure of t2g-orbitals, which are split by the crystal field, is rather similar in YTiOa, 
GdTiOs, and SmTiOa: apart from some numerical differences, the coefficients of the 
expansion over the basis function have similar values and similar phases. The same 
tendencies are clearly seen in the RT structure of YTiOs after the interchanging of the 
middle and highest t2g-levels. Nevertheless, the structure of the t2g-levels in LaTiOa 
appears to be different. 

The behavior of transfer integrals in the NN bonds is explained in table O 
The transfer integrals are presented in the local coordinate frame of orbitals, which 

Table 2. Transfer integrals in the bonds 1-2 and 1-3 (measured in meV) in the 
local coordinate frame corresponding to the diagonal representation of the crystal-field 
splitting. The positions of the atomic sites are explained in Fig. [TJ Transfer integrals 
in other bonds can be obtained from ti2 and ii3 using the symmetry operations of the 
space group D^f^. 

compound ti2 ti3 

/ -8 91 13 \ / -29 71 -12 \ 
YTiOs (RT) -4 -66 71 82 -38 

V 193 11 -32 / V -12 -38 102 / 

/ -36 26 96 \ / -49 -12 64 \ 
YTiOs (LT) 185 -8 39 { ^7 -33 

V -26 -7 -68 / V 64 -33 114 / 

/ -48 54 85 \ / -34 13 80 \ 
GdTiOa 171 6 23 13 98 -24 

V -56 -20 -79 / V 80 -24 143 / 

/ -83 43 81 \ / -10 1 101 \ 
SmTiOa 159 39 27 1 114 -12 

V -45 -23 -84 / V 101 -12 143 / 

/ -62 130 94 \ / 142 -70 78 \ 
LaTiOg -14 -5 94 -70 90 53 

V 76 76 -69 / V 78 53 53 / 



diagonalize the crystal field (see tabled]). 

The screened on-site Coulomb interactions for the t2g-band are expressed through 
the 3x3x3x3 matrices, in the basis of three Wannier orbitals [8[ HOl Ej. For the 
distorted compounds, the structure of these matrices is rather complex and may involve 
many independent parameters. These matrices are directly used in the next section for 
the analysis of SE interactions without any additional approximations or simplifications. 
Nevertheless, just for explanatory purposes in this section, the full matrix Ht/a/j-y^H was 
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fitted in terms of two Kanamori parameters [13]: tfie Coulomb interaction U between tlie 
same t2g-orbitaIs and tlie excliange interaction J'. The fitting implies that Ht/a/j^^H has 
the same symmetry as in the spherical environment of isolated atoms. For example, in 
the process of fitting, it was assumed that the Coulomb interaction W between different 
t2g-orbitals is related to U and J' by the identity W = IA — 2 J [Mj, and the parameter U 
is the same for all three t2p-orbitals. The results of such fitting are shown in table [31 The 

Table 3. Results of fitting of the effective Coulomb interactions in terms of two 
Kanamori parameters: the intraorbital Coulomb interaction hi and the exchange 
interaction . All energies are measured in eV. 



compound 


U 


J 


YTiOa (RT) 


3.45 


0.62 


YTiOg (LT) 


3.40 


0.62 


GdTiOs 


3.28 


0.62 


SmTiOa 


3.29 


0.62 


LaTiOa 


3.20 


0.61 



screening of the Coulomb interaction lA is sensitive to the local environment in solids: 
generally, U is larger for the more distorted YTiOa and smaller for the least distorted 
LaTiOa [8]. On the contrary, the exchange interaction J is less sensitive to the details 
of the screening and close to the atomic limit [T5] . 

4. Superexchange Interactions in the One-Electron Approximation 

The superexchange interaction in the bond i-j is related to the gain of the kinetic energy, 
which is acquired by an electron residing at the atomic site i in the process of virtual 
hoppings in the subspace of unoccupied orbitals at the atomic site j, and vice versa 
[ini [T7j. Let us consider first the approximation adopted in [TOl [18], where the energy 
gain caused by virtual hoppings in the bond i-j was computed in the following way: 

Namely, we start with the lattice of isolated atoms, each of which accommodates one 
electron, and describe the ground-state wavefunction of such a reference system by the 
single Slater determinant G. Since in the atomic limit there is only one t2g-electron 
per one Ti-site, this is essentially one-electron problem and all many-electron effects 
emerge only in the process of virtual hoppings. Typically, this justifies the use of the 
single-determinant approximation for G in the case of titanates By denoting the 
occupied one-electron orbitals at the sites i and j as ipi and ipj, respectively, the Slater 
determinant G for the bond i-j takes the following form: 

|G(1, 2)) = 4^ b.(l)^,(2) - ip^{2)ip,{l)] , 



tij 



E 

M 



^\jM){jM\lP, 



'jM 



G 



(4) 
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where the symbols "1" and "2" stand for the coordinates of two electrons, associated 
with the sites i and j. After that we treat the transfer integrals tij as a perturbation. 
Since the Coulomb repulsion U (and W) is large, it is sufficient to consider the excited 
configurations accommodating only two t2g-electrons, which contribute to the second 
order perturbation theory with respect to tij. The eigenvalues and eigenvectors of such 
excited configurations at the site j are denotes as EjM and \jM), respectively. They 
are obtained after the diagonalization of the Coulomb interactions, 

in the basis of all possible two-electron Slater determinants at the site j . For six t2g spin- 
orbitals, there are 6(6— 1)/2=15 such determinants [lO]. Thus, each |jM) is constructed 
from several Slater determinants and takes into account the correct multiplet structure 
in the atomic limit. Examples of the atomic multiplet structures corresponding to the 
exciting two-electron configurations can be found in [10]. in (jl]) is a projector 
operator, which enforces the Pauli principle and suppresses any transfer of an electron 
into the subspace of occupied orbitals at the site j. In practice, projects |jM) into 
the subspace spanned by the Slater determinants of the form 

|G,^,(1,2)) = i=[^,(l)^,-(2)-^,(2)^,-(l)], 

where ipj is the occupied orbital at the site j and ipj is any orbital residing at the site 
J, which can be reached from ipi by the transfer integrals tij. 

Thus, the expression for the energy gain combines elements of the one-electron 
approximation for \G) with the exact many-electron treatment for the excited two- 
electron states. On the one hand, since \jM) is a combination of several Slater 
determinants, the method takes into account some many-electron effects in the atomic 
limit, which may have interesting consequences on the magnetic properties of i^TiOsJfl 
On the other hand, the form of \G) is restricted by the single Slater determinant. In 
this sense, this is an one-electron approach. That is why, in the following we will refer 
to the expression (jl]) as to an one-electron approximation for the SE interactions. The 
multi-determinant analog of (jl]) will be considered in Sec. [5l 



4--1- Crystal- Field Theory for the Superexchange Interactions at T=0 

In the crystal-field theory, it is assumed that the form of the occupied orbitals {v^i} is 
totally controlled by the crystal field, which is much larger than the energy gain caused 
by the virtual hoppings (jl]) as well as the energies of thermal fiuctuations. Then, in 

^ For example, in the case of the AFM spin ahgnment in the bond i-j, the excited configuration "fl 
is decomposed into two-electron singlet and triplet states, that leads to the additional energy gain in 
equation (j4]). This effect additionally stabilizes the AFM interactions, as is manifested for example in 
somewhat stronger canting of spin magnetic moments away from the collinear FM alignment in the 
ground state of YTiOs in comparison with results of the mean- field Hartree-Fock approximation [SlllOj. 
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the one-electron approximation, the isotropic spin couphngs Jij is related to the energy 
gains (jl]) by the following expression 

2 7- • = T^'^ - T^.^ 

where T^j^ = T^ip], ipj) and T^^^ = ^(v^I; V^j)? the indices t and | explicitly show the 
spin states of the occupied orbitals y?, and ipj. The values of J^j between the nearest and 
some next nearest neighbors are listed in table HI It also shows the corresponding type of 
the magnetic ground state expected in the i^TiOs compounds. The magnetic transition 



Table 4. Results of the crystal-field theory for RTiO^: isotropic superexchange 
interactions {Jij, measured in meV), corresponding magnetic transition temperatures 
(7c,N, measured in K) obtained in the random phase and mean-field approximation 
(results of the mean-field approximation are shown parenthesis), and the type of the 
magnetic the ground state (GS, where "F" stands for the ferromagnetic state and "A" - 
for the A- type antiferromagnetic state). The positions of the atomic sites are explained 
in Fig. [T] Depending on the magnetic ground state, the notations Tq and Tn stand 
for the Curie and Neel temperature, respectively. 



compound 


Jl2 


Jl3 


J23 


■723' 




GS 


YTiOs (RT) 


3.23 


0.45 


0.03 


-0.16 


74 (154) 


F 


YTiOs (LT) 


2.93 


-0.14 


0.03 


-0.16 


85 (144) 


A 


GdTiOa 


2.79 


0.68 


0.07 


-0.17 


79 (138) 


F 


SmTiOa 


0.68 


1.60 


0.07 


-0.03 


41 (68) 


F 


LaTiOa 


1.41 


-4.88 


0.30 


0.13 


93 (164) 


A 



temperature is estimated in the random phase approximation (RPA, Appendix A ). 

Thus, in the CF theory, the ground state appears to be FM in the case of GdTiOs 
and SmTiOa, and A- type AFM in the case of LaTiOs. The ground state of YTiOa 
depends on the crystal structure: FM and A-type AFM for the RT and LT structure, 
respectively. Nevertheless, as we will see in the next section, both structures yield the 
same type of the noncoUinear magnetic ground state, which combines elements of the 
FM and A-type AFM ordering within one magnetic structure, when the relativistic 
SO interaction is taken into account. To certain extent the FM interactions in i^TiOa 
can be also stabilized by considering explicitly the e^-band of these compounds [5J. 
Nevertheless, at the present stage the situation is not completely clear because the 
same effects would act against the G-type AFM ground state in the case of SmTiOa 
and LaTiOs. 

The magnetic interactions in YTiOa (RT) are in reasonable agreement with the 
ones obtained in the Hartree-Fock approximation for the Hubbard model ([T]) with the 
same parameters of the crystal structure ( Ji2= 3.2 ^ 3.9 meV and Ji3= 1.0 ^ 1.2 meV, 
depending on the magnetic state [lO])- To certain extend, the same is true for LaTiOa, 
although since the CF splitting is small and allows for the additional change of the 
orbital ordering in each magnetic state [17], the parameters of interatomic magnetic 
interactions obtained in the Hartree-Fock approximation for LaTiOs exhibit a strong 
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dependence on the magnetic state in which they are calculated (Ji2= 1.0 -i- 4.5 meV 
and Ji3= —1.2 -. — 4.9 meV [10]). Nevertheless, results of the CF theory (table H]) are 
within the parameters range obtained previously in the Hartree-Fock approximation for 
LaTiOs [TO]. 

RPA considerably reduces the magnetic transition temperature in comparison 
with the mean-field approach. This effect is particularly strong in YTiOa due to 
the quasi-two-dimensional character of the SE interactions, which according to the 
Mermin- Wagner theorem suppress the long-range magnetic order at finite T [19jFI 




Nevertheless, J13 is finite and the reduction of the magnetic transition temperature 
has only logarithmic dependence of the anisotropy g = IJ13/J12I of NN interactions: 
^c,N ~ ■~(liifi')~^ [20] El Therefore, the transition temperature remains finite. Moreover, 
the Curie temperature is overestimated by factor two even in the case of YTiOs. The 
correct G-type AFM ground state is reproduced neither in SmTiOs nor in LaTiOs. 

4-2. Relativistic Spin- Orbit Interaction and the Effective Spin Hamiltonian 

In the relativistic generalization of the crystal-field theory, the occupied orbital ipi at each 
Ti-site is obtained after the diagonalization of the one-electron Hamiltonian hf^ + hf^, 
which combines the crystal field and the relativistic SO interaction hf^ = (^j/2)(crj, Ij), 
where C,i ~ 20 meV is related to the spherical part of the one-electron potential [8]. The 
Hamiltonian is constructed in the basis of six t2g spin-orbitals. The lowest eigenstate 
obtained after the diagonalization is the Kramers doublet, whose eigenvectors can be 
formally denotes as if] and iff. Then, we consider a liner combination 



and find the coefficients c} and of from the condition that the averaged spin moment, 
ef°'={ipf°'\a-i\ipf°') corresponding to ipf"', has the maximal projection along the ±a = 
±x, zby, and ±z axes in the orthorhombic coordinate frame^l Then, we use these orbitals 
in the expression (jlj) for the energy gain, and calculate 36 parameters T{Lpf°', Lp^^) 
corresponding to all possible combinations of ±a and ±6 at the sites i and j. In 
the one-electron (mean-field) approximation, the magnetic part of T{Lpf"', '■Pf'') should 

^ Strictly speaking, this behavior does not seem to be consistent with the experimental inelastic 
neutron scattering data, which are typically interpreted in terms of the three-dimensional isotropic 
Heisenberg model [211. However, at present there is no clear consensus on this matter, neither on 
theoretical nor on experimental side. For example, the experimental orbital ordering pattern determined 
in [22] is more consistent with the anisotropic (quasi-two-dimensional) structure of interatomic magnetic 
interactions |101l23j. Thus, the problem requires additional study, both on theoretical and experimental 
sides. 

^ Moreover, the dependence Tc.n ^ — (lng)~^ can be further modified by the longer range interactions 
J23 and J23' between neighboring ab-planes. 

^ The actual procedure was based on the numerical maximization of the function 



(e±^e±'^)/(e±^e±°)l/2, where = (±1,0,0), = (0,±1,0), and = (0,0, ±1). The proce- 





i 1 
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correspond to the energies, which are obtained from ([2]) by replacing the Pauh matrices 
&i by the unit vectors e^, describing different directions of spin at the site i. In 
total, the magnetic part of T{(pf°', cp^^) is characterized by 9 independent parameters, 
which constitute the tensor Aij of interatomic magnetic interactions. The values of 
Aij for two inequivalent NN bonds are given in table [5] and the decomposition of 
Aij into Jij, dij, and fij is presented in table [61 One can clearly see that the 

Table 5. Tensors of magnetic interactions associated with the nearest- neighbor 
bonds 1-2 and 1-3 (measured in nieV) as obtained in the crystal-field theory with 
the relativistic spin-orbit interaction. The positions of the atomic sites are explained 
in Fig. [TJ The magnetic interactions in other nearest-neighbor bonds can be obtained 
from Ai2 and A13 using the symmetry operations of the space group D^f^ . 



compound A12 j4i3 



YTiOs (RT) 


/ -3.11 
0.22 
V 0.07 


0.22 
-3.11 
-0.12 


0.06 \ 
0.30 
-3.26 ) 


/ -0.38 
-0.07 
V 0.25 


-0.07 
-0.39 
0.37 


-0.25 \ 
-0.37 
-0.46 / 


YTiOs (LT) 


/ -2.89 
-0.01 
V 0.41 


0.26 
-2.86 
-0.37 


-0.33 \ 

0.48 
-2.97 ) 


/ 0.19 
-0.01 
V 0.10 


-0.01 
0.18 
0.31 


-0.10 \ 
-0.31 
0.17 / 


GdTiOa 


/ -2.79 
-0.09 
V 0.70 


0.40 
-2.66 
-0.42 


-0.59 \ 

0.62 
-2.84 ) 


/ -0.62 
-0.03 
V 0.22 


-0.03 
-0.62 
0.34 


-0.22 \ 
-0.34 
-0.65 / 


SmTiOa 


/ -0.66 
-0.41 
V 1.22 


0.72 

-0.56 
-0.72 


-1.07 \ 
0.93 

-0.72 ) 


/ -1.49 
-0.05 
V 0.13 


-0.05 
-1.57 
0.34 


-0.13 \ 
-0.34 
-1.58 / 


LaTiOa 


/ -0.86 
1.84 
V 0.24 


-1.22 
-1.20 
-1.71 


-1.24 \ 

0.82 
-0.93 ) 


/ 4.81 
0.34 
V 0.56 


0.34 
4.21 
-3.36 


-0.56 \ 
3.36 
4.43 / 



magnetic interactions d^j and of the relativistic origin steadily increase in the direction 
YTiOs^GdTiOs^SmTiOa^LaTiOs, which is quite consistent with the decrease of the 
CF splitting in the same direction (see tabled]). Nevertheless, it is interesting to note 
that the relativistic SO interaction also contributes to the isotropic interactions Jij, and 
all of them decrease (with some exception of LaTiOs) after taking into account the SO 
interaction in comparison with results of the pure CF theory in table HI 

The directions of the magnetic moments in the ground state and the values of 
the magnetic transition temperature in the mean- field approximation ( Appendix B ) are 
summarized in table [3 As expected, the magnetic ground state is noncoUinear. The type 
of the magnetic ground state for the space group D^l can be formally denoted as X-Y-Z, 
where X, Y, and Z is the magnetic structure (F, A, C, or G) formed by the projections 
of the magnetic moments onto the orthorhombic axes a, b, and c, respectively. By 
comparing the values of the magnetic transition temperature with the ones reported in 
table m one can clearly see that the SO interaction tends to additionally increase Tq, 
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Table 6. Isotropic Heisenberg interactions (Jy ), antisymmetric Dzyaloshinsky-Moriya 
interactions (dij), and symmetric anisotropic (fy) interactions associated with the 
nearest-neighbor bonds 1-2 and 1-3 as obtained in the crystal-field theory with the 
relativistic spin-orbit interaction. All interactions are measured in meV. The positions 
of the atomic sites are explained in Fig. [T] 



compound 


Jl2 


di2 




Tl2 




YTiOg (RT) 


3.16 


( "^'^^ ^ 
-0.01 

V ) 


/ 0.05 
0.22 
V 0.07 


0.22 
0.05 
0.10 


0.07 \ 
0.10 
-0.10 / 


YTiOs (LT) 


2.90 


/ -0.42 \ 

-0.37 
V -0.13 j 


/ 0.02 
0.12 
V 0.04 


0.12 
0.05 
0.06 


0.04 \ 
0.06 
-0.07 / 


GdTiOa 


2.76 


1 -0.52 ^^ 

-0.65 
V -0.25 j 


/ -0.02 
0.15 
V 0.05 


0.15 
0.10 
0.10 


0.05 \ 
0.10 
-0.08 / 


SmTiOa 


0.65 


/ -0.83 ^^ 

-1.15 
V -0.56 J 


/ -0.02 
0.16 
V 0.07 


0.16 
0.09 
0.10 


0.07 \ 
0.10 
-0.07 / 


LaTiOs 


0.99 


( ""^'^^ ^ 
-0.74 

V 1.53 ) 


/ 0.14 
0.31 
V -0.50 


0.31 
-0.20 
-0.45 


-0.50 \ 
-0.45 
0.06 / 




compound 












YTiOs (RT) 


0.41 


/ 0.37 \ 

-0.25 
V ) 


/ -0.03 
-0.07 
\ 


-0.07 
-0.02 



° ) 

0.05 / 


YTiOs (LT) 


-0.18 


1 0.31 \ 
-0.10 

V ) 


/ 0.01 
-0.01 
V 


-0.01 




" ) 

-0.01 / 


GdTiOa 


0.63 


/ 0.34 ^^ 
-0.22 

V J 


/ 0.01 
-0.03 
\ 


-0.03 
0.01 



» ) 

-0.02 / 


SmTiOa 


1.55 


/ 0.34 ^^ 

-0.13 
V J 


/ 0.06 
-0.05 
\ 


-0.05 
-0.03 



" ) 

-0.03 / 


LaTiOa 


-4.48 


1 -3.36 \ 
-0.56 

V ) 


/ -0.33 
-0.34 
\ 


-0.34 
0.27 



» ) 

0.06 / 



despite the fact that all isotropic interactions Jij decrease. Moreover, the SO interaction 
alone does not solve the problem of the magnetic ground state of LaTiOa, which is 
expected to be of the C-F-A type and does not includes the experimentally observed 
G-component. The magnetic ground state of other compounds is of the G-A-F type. 
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Table 7. Magnetic ground state, direction of magnetization in the ground state 
and the magnetic transition temperature (Tc,n, measured in K) in the mean- field 
approximation as obtained in the crystal-field theory with the relativistic spin-orbit 
interaction. The vector of magnetization is referred to the site 1 in Fig. [TJ Similar 
vectors at other Ti-sites of the primitive cell are obtained by applying the symmetry 
operations of the space group D^^. 



compound 


ground state 


direction of magnetization 


7c, N 


YTiOa (RT) 


G-A-F 


(-0.02,-0.54, 0.84) 


164 


YTiOa (LT) 


G-A-F 


(-0.05,-0.89, 0.45) 


148 


GdTiOa 


G-A-F 


(-0.12,-0.38, 0.92) 


150 


SmTiOg 


G-A-F 


(-0.35,-0.20, 0.91) 


94 


LaTiOa 


C-F-A 


(-0.21,-0.47, 0.86) 


197 



Particularly, in YTiOa, there is a considerable weight of the both A- and F-components 
parallel to the orthorhombic b- and c-axes, respectively. The weight of the G-component 
parallel to the a- axis is negligibly small. The main difference between the RT and LT 
structures is in the relative weight of the A- and F-components (the formed is larger in 
the LT structure). On the other hand, in SmTiOa there is a substantial weight of the 
G-type AFM component along the a-axis. The magnetic structures obtained for YTiOa 
(RT) and LaTiOs are consistent results of the Hartree-Fock calculations, which were 
considered in the previous publication [10]. 

4-3. Thermal Fluctuations of the Orbital Degrees of Freedoms 

In this section we consider the effect of thermal fluctuations of the orbital degrees of 
freedom on interatomic magnetic interactions between the spins. For these purposes, 
it is convenient to work in the local coordinate frame corresponding to the diagonal 
representation of the CF Hamiltonian hf^ at each atomic site. Then, we assume that 
for each projection of spin, the three-component occupied t2g-orbital at the site i can 
be presented in the form of the real vector 

^ cos 6i ^ 

Vi = sin 9i cos 0j , (5) 
^ sin 9i sin (pi j 

in the basis of three CF orbitals listed in tabled], where < 6i < tt /2 and < (pi < n 
due to the invariance of both crystal field and pair interactions (jl]) with respect 
to the inversion Vi — Wj. In these notations, the point 6i=(pi=0 corresponds to 
the lowest CF orbital. This situation was already considered in Sec. 14.11 in the 
limit of large CF splitting. However, since the CF splitting is finite, other orbital 
configurations can contribute to the thermodynamic averages of physical quantities at 
elevated temperatures. Thus, similar to Sec. 14.11 we assume that the CF splitting 
is much larger that the energy gain (jl]) caused by the virtual hoppings, but can 
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become comparable with the energies of thermal fluctuations, and consider the flnite- 
temperature extension of the CF theory. Then, the thermal average of (jl]) is given 
by 

%,{T) = j j dQ,w{e„(i)j)T{e,,<i)„e„<Pj), (6) 

where 



27r i>n/2 

dVli = — I d(f)i / sin 9id6i, 
Jo 



1 

2^ 



Zi = dVti exp 



and T{6i, (pi, 6j, is the pair interaction (jlj) constructed from the occupied t2g-oi'bitals 
of the form ^ for either ferromagnetic (|t) or antiferromagnetic {][) conflgurations of 
spins (in the following denoted as 7^j^ and 7^]^, respectively). The numerical integration 
in ([6]) was performed by using the the Metropolis algorithm [2H [25]. An example of 
the temperature dependencies of 7^j^(T) and T^-\T) for YTi03 and LaTiOa is shown in 
Fig. [2]0 Then, using 7^j^(T) and T^-^(T), one can evaluate the temperature dependence 
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Figure 2. Thermal averages of pair interactions Tij{T) for the ferromagnetic (tt) 
and antiferromagnetic (ti) configurations of spins in the bonds 1-2 and 1-3 calculated 
for the low-temperature structure of YTiOa (left) and LaTiOs (right). The atomic 
positions are explained in Fig. [1] 

^ The behavior of GdTiOs and SmTiOa is rather similar to that of YTiOa and not discussed here. 
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of interatomic interactions between the spins, 2Jjj(T)=7^j^(T)— 7^j^(T), averaged over 
all orbital configurations, and self-consistently solve the equation for the magnetic 



transition temperature in RPA ( Appendix A ), where the value of Tc,n is used as the 
argument of Jij(T), and the procedure is repeated until reaching the self-consistency 
with respect to Tc,n- Then, we obtain the following values of the Curie temperature 
Tc= 60, 68, and 37 K for YTiOg (RT), GdTiOg, and SmTiOg, respectively. YTiOs (LT) 
and LaTiOs are expected to develop the A-type AFM order with the Neel temperature 
Tn= 64 and 52 K, respectively. Thus, the thermal fluctuation systematically decreases 
the values of the magnetic transition temperature (by 10-44 %). As expected, the largest 
change is observed in LaTiOs, which has the smallest CF-splitting, and in YTiOa, due to 
the quasi-two-dimensional character of interatomic interactions. Nevertheless, SmTiOs 
remains FM and LaTiOs ^ A-type AFM, contrary to the experimental data. 



5. Multi-Determinant Approach for Superexchange Interactions 

By summarizing results of the previous sections, we note the following shortcomings of 
the one-electron approach: 

• In the case of LaTiOa and SmTiOa, it fails to predict the correct G-type AFM 
ground state. Although small G-type AFM component (along the a-axis) is 
expected in the theoretical magnetic ground state of SmTiOa after including the 
relativistic SO interaction, that of LaTiOs does not involve the G-type AFM 
arrangement. 

• Even for the FM compounds YTiOs and GdTiOs, the magnetic transition 
temperature is typically overestimated by factor two. To certain extent, it can 
be reduced by thermal fluctuations of the orbital degrees of freedom. On the other 
hand, the relativistic SO interaction acts in the opposite direction and additionally 
increases Tq (at least on the level of mean- field approximation). 

In this section, we investigate whether these problems can be resolved by considering 
the many-electron effects. 

As was pointed out in the beginning of Sec. HI the main drawback of all previous 
considerations was that, although the expression (jlj) for the energy gain takes into 
account the many-electron effects in the intermediate configurations t^g^ which can 
be reached in the process of virtual hoppings, it was combined with the one-electron 
approximation for the ground-state wavefunction G in the atomic limit. The purpose 
of this section is to go beyond this one-electron approximation and to clarify the role 
played by the many-electron effects in the problem of SE interactions. Namely, for each 
bond i-j, we construct the complete basis of two-electron Slater determinants of the 
form, 

|S(1,2)) = i= [^,(1)^,(2) -^^.(2)^,(1)] , (7) 

where ipi denotes the one-electron spin-orbital residing at the site i. In practice, these 
orbitals were obtained from the diagonalization of the CF Hamiltonian hf^ (tabled]). In 
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tji + 



S' (8) 



total, there are 6 such spin-orbitals and 36 Slater determinants of the form ([7]). Then, 
for each bond, we calculate matrix elements of the pair interactions 

#,|jM)(jM|# / 

EjM 

in the basis of these Slater determinants, and combine them with the matrix elements 
of one-electron operators of the crystal field hf^ + h^^ and (optionally) the SO 
interaction /i?^ + h^'^ . The corresponding 36x36 Hamiltonian matrix is denoted as 
Hij=\^Hfj^' Then, the spectrum of two-electron states in each bond is obtained 
after the diagonalization of Hij. Thus, the basic idea (and approximation) behind 
this treatment is that the entire lattice is divided into "molecules", and the electronic 
structure of each "molecule" can be considered independently from other "molecules". 
We will mainly focus on the origin of AFM correlations in LaTiOs. 

First, let us discuss results without relativistic SO interaction. The spectrum of 
two-electron states, obtained after the diagonalization of Hij for the bonds 1-2 and 1-3 
is shown in Fig. [3l There is a drastic difference of LaTiOs from other compounds. In 
YTiOs, SmTiOs, and GdTiOs, four low-energy levels, including one spin-singlet and 
three degenerate spin-triplet states, are clearly separated from the next states by an 
energy gap of at least 100 meV, corresponding to the value of the CF splitting (tabled]). 
However, in LaTiOa, the singlet-triplet splitting in the low-energy part of the spectrum is 
at least comparable with energy gap, separating these levels from the next two-electron 
states. For example, in the bond 1-3, the singlet-triplet splitting is about 20 meV. 
However, the next level is located only 4 meV higher than the triplet state. Thus, the 
crystal-field theory, although applicable for YTiOs and high-temperature structures of 
SmTiOs and GdTiOs, definitely breaks down in the case of LaTiOa, where the splitting 
of the two-electron levels caused by the SE effects in the low-energy part of the spectrum 
is at least comparable with the CF splitting. 

The parameter of the isotropic exchange coupling in the bond i-j can be expressed 
through the energies of the low-lying spin-singlet {Es) and spin-triplet {Et) states, 

Jij = ^ {Es - Et) , 

where the prefactor 1/4 stands for 5^, which according to ([3]) is already included to 
the definition of Jij. The obtained parameters are listed in table [H] together with the 
values of the magnetic transition temperature Tc,n estimated in the random phase 
approximation ( Appendix A ). By comparing them with the values obtained in the 
one-electron approximation (table Hj), one can clearly see that the two-electron effects 
tend to additionally stabilize the FM interactions. Generally, the values of all FM 
interaction Jij increase in the two-electron approach, while the magnitude of the AFM 
interactions J13 in LaTiOs and in the LT phase of YTiOs decreases. Particularly, Jis 
nearly vanishes in the case of YTiOs (LT), resulting in the two-dimensional character 
of NN interactions. However, due to the next NN interactions between the planes, 
Tn remains finite. Nevertheless, only in YTiOs (LT) the two-electron effects tend to 
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Figure 3. Spectrum of two-electron states obtained after the diagonalization of pair 
and crystal-field interactions in the complete basis of Slater determinants separately 
for the bonds 1-2 and 1-3 (the atomic positions are explained in Fig.[T]). 



decrease the magnetic transition temperature. This behavior is solely related to the 
quasi-two-dimensional character of interatomic magnetic interactions. 

Thus, it seems unlikely that the two-electron effects alone will stabilize the G-type 
AFM ground state in the case of LaTiOa. Therefore, we investigate the last possibility 
related to the relativistic SO interaction. For the most of the considered systems, the SO 
interaction is small in comparison with the CP splitting and does not significantly change 
the distribution of the two-electron states. The typical example for the LT structure 
of YTiOa is shown in Fig. HI which is practically identical to the spectrum without the 
SO interaction (Fig. [3]). For the FM systems, the splitting of the low-lying spin-triplet 
states by the SO interaction is very small. Nevertheless, below we will see that this 
splitting may have rather interesting consequences on the magnetic properties. Again, 
one clear exception is LaTiOs, where the strength of the SO interaction = 21 meV) is 
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Table 8. Parameters of isotropic exchange coupling {Jij, measured in meV) derived 
from the singlet-triplet splitting of the low-energy levels, which were obtained from 
the diagonalization of pair and crystal-field interactions in the complete basis of two- 
electron Slater determinants separately for each bond of the system; corresponding 
magnetic transition temperatures (Tcn, measured in K) obtained in the random phase 
approximation; and the type of the magnetic ground state. Depending on the magnetic 
ground state, the notations Tq and Tn stand for the Curie and Neel temperature, 
respectively. 



compound 


Jl2 
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J23 
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0.98 


0.07 
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Figure 4. Spectrum of two-electron states obtained for the low-temperature structure 
YTiOa (left) and LaTiOa (right) after the diagonalization of pair, crystal-field 
and relativistic spin-orbit interactions in the complete basis of Slater determinants, 
separately for the bonds 1-2 and 1-3 (the atomic positions are explained in Fig.[T]). 

comparable with the CF sphtting and the energy gain caused by the virtual hoppings, 
and only in LaTiOs the spectrum of two-electron states changes significantly when the 
relativistic SO interaction is taken into account. In order to understand the character 
of magnetic interactions, we calculate the parameters of interatomic spin correlations in 
the ground state 

where \1/ is the two-electron wavefunction for the bond i-j, corresponding to the lowest 
eigenvalue of Hij. These parameters are listed in table [9l The diagonal matrix elements 
{a=b) describe the longitudinal spin correlations. If ((^"(j") > 0, the correlations are FM. 
If {d'°'a'j) < 0, they are AFM. The off-diagonal elements {a^h) describe the transverse 



Superexchange Interactions in Orthorhombically Distorted Titanates RTiOs (R= Y, Gd, Sm, and La)19 



Table 9. Parameters of interatomic spin correlations in the ground state calculated 
separately for each bond of the system after the diagonalization of the two-electron 
Hamiltonian including pair, crystal field, and relativistic spin-orbit interactions. The 
atomic positions are explained in Fig. [1] 



compound 








1 






YTiOs (RT) 


/ 0.03 
-0.92 
V -0.31 


-0.91 
0.08 
-0.35 


-0.35 \ 
-0.29 
0.83 ) 


/ -0.40 
0.83 
V 0.29 


0.83 
0.21 
0.44 


-0.29 \ 
-0.44 
0.77 / 


YTiOa (LT) 


/ 0.10 
-0.89 
\ -0.42 


-0.95 


-0.23 


-0.22 \ 
-0.43 
0.85 / 


/ -0.91 
0.25 
\ 0.28 


0.25 
-0.14 
0.93 


-0.28 \ 
-0.93 
-0.07 / 


GdTiOa 


/ 0.30 
-0.91 
V -0.16 


-0.80 
-0.19 
-0.53 


-0.47 \ 
-0.29 
0.79 ) 


/ -0.41 
0.86 
V 0.21 


0.86 
0.33 
0.33 


-0.21 \ 
-0.33 
0.88 / 


SmTiOa 


/ 0.13 
-0.92 
V 0.12 


-0.51 
-0.23 
-0.78 


-0.79 \ 
-0.06 
0.48 ) 


/ -0.66 
0.71 
V 0.06 


0.71 
0.65 
0.14 


-0.06 \ 
-0.14 
0.95 / 


LaTiOg 


/ -0.23 
-0.44 
V 0.71 


-0.67 
0.28 
0.08 


0.42 \ 
0.59 
0.27 ) 


/ -0.62 
-0.01 
\ 


-0.01 
-0.55 
-0.16 


0.16 
-0.90 / 



correlations, which are related to the noncoUinear spin arrangement in the bond i-j. 
Generally, the longitudinal correlations are strongly anisotropic. For example, within 
one bond the correlations can easily become either FM or AFM, depending on the 
direction in the orthorhombic lattice. The transverse correlations, which are solely 
caused by the relativistic SO interaction, can be also strong and comparable with the 
longitudinal ones. 

For the FM systems, this picture may substantially differ from the one without the 
SO interaction. However, it should be noted that the details of the ground state in the 
FM case are determined by the small splitting of the low-lying spin-triplet levels by the 
SO interaction. For the most of the systems this splitting is small and typically varies 
from few hundredths of meV to one meV. For example, for the FM bond 1-2 in YTiOa 
(LT), the second and third levels are split from the lowest one by 0.56 meV and 0.59 
meV, respectively. Therefore, when the temperature exceeds T ~ 0.6 meV/fc^ ~ 7 K, 
the proper picture for the magnetic interactions in this compound would correspond to 
the (thermal) average || ((3""(3"2) || over the low-lying "spin-triplet" states. For example, 
for the 1-2 bond in YTiOs (LT) this procedure yields the following tensor 

0.03 \ 
0.04 . 

0.32 y 
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Similar results can be obtained for other systems forming the FM bonds. After the 
averaging, the longitudinal spin correlations become isotropic, while all transversal 
correlations are considerably smaller. This behavior is consistent with the form of 
the interatomic magnetic interactions tensor Aij obtained on the level of one-electron 
approximation (table [5]). Another effect, which can mix the two-electron states within 
the lowest manifold of certain bond is related to the virtual hopping in neighboring 
bonds, which involve one of the atomic sites of the original bond. Note that the many- 
electron effects will mix the lowest CF orbital at certain atomic site with other states. 
The symmetry of these states as well as the magnitude of the mixing will be generally 
different for different bonds. Since each atom participate in several bonds, this effect 
will lead to the addition mixing of the low- lying (spin-triplet) states obtained from the 
diagonalization of the two-electron Hamiltonian separately for each bond of the system. 
Nevertheless, the quantitative estimate of this effect requires a more rigorous solution 
of the many-electron problem, which is beyond the scopes of the present work. 

In the case of LaTiOs, all longitudinal correlations in the bond 1-3 are 
antiferromagnetic, while the transversal correlations are small. The correlations in the 
bond 1-2 are strongly anisotropic: if yy- and 2;z-components favor the FM coupling, 
the xx-correlations are antiferromagnetic. The splitting between the lowest and the 
next two two-electron levels in the bond 1-2 is about 4 meV, which is considerably 
larger than the splitting obtained in other (more distorted) compounds. Therefore, 
the effect is expected to be more robust against the thermal fluctuations mixing the 
spin-triplet states. Strong transversal correlations are also expected. Thus, the xx- 
correlations appears to be antiferromagnetic simultaneously in the bonds 1-2 and 1-3, 
being consistent with the G-type AFM structure. On the basis of this analysis, we 
expect that if the G-type AFM order took place in LaTiOa, it would be more likely 
developed by the x- (a-) projections of the magnetic moments in the orthorhombic 
coordinate frame. 

6. Summary and Conclusions 

Starting from the multiorbital Hubbard Hamiltonian for the t2g-bands of orthorhom- 
bically distorted titanates i^TiOs {R= Y, Gd, Sm, and La), where all the parameters 
were derived from the first-principles electronic structure calculations, we considered 
different levels of approximations for the construction of the spin-only superexchange 
model. Namely, after considering the conventional crystal-field theory, where all inter- 
atomic magnetic interactions are solely determined by occupied t2g-orbitals, that are 
split off by the crystal distortion, we consecutively incorporate the effects of the rel- 
ativistic spin-orbit interaction, thermal fluctuations of the orbital degrees of freedom, 
and many-electron effects related to the virtual electron hoppings in the bonds. 

Even in the conventional CF theory, the interatomic magnetic interactions appear 
to be extremely sensitive to the details of the orbital structure, and small change of the 
orbital structure, caused by either crystal distortions or thermal fluctuations, can have 
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a profound effect on the magnetic properties. 

Particularly, the use of the room- and low-temperature structure for YTiOa provides 
rather different sets of parameters of the spin Hamiltonian. The additional distortion 
in the low-temperature structure tends to weaken the ferromagnetic interactions: it 
substantially reduces the FM coupling J12 in the orthorhombic ab-plane and makes 
the coupling along the c-axis weakly AFM. Similar tendency was reported for LaTiOs, 
where the use of the low-temperature structure [13] could help to stabilize the G-type 
AFM state [5l [71 [26] . Apparently, some of the problems encountered in the present 
work, such as the incorrect magnetic ground state in SmTiOa and the overestimation 
of the Curie temperature in GdTiOs, may be resolved by using the low-temperature 
structural data for these compounds, which are not available today. 

Moreover, the thermal fluctuations of the of the orbital degrees of freedom near 
the CF configuration can substantially reduce the value of the magnetic transition 
temperature (up to 40%). 

The relativistic spin-orbit interaction, which is responsible for the appearance 
of anisotropic and antisymmetric Dzyaloshinsky-Morita interactions, leads to the 
noncollinear magnetic alignment. The role of the relativistic effects increases in the 
direction YTiOs— i>GdTi03— s^SmTiOs— s^LaTiOs, when the crystal distortion decreases. 

The crystal-field theory, although applicable for YTiOa and high-temperature 
structures of GdTiOs and SmTiOs, definitely breaks down in the case of LaTiOa, 
which has the smallest CF-splitting and where other factors, such as the relativist SO 
interaction and many-electron effects in the bonds, start to play an important role. 
Particularly, we found that the combination of the latter two factors could explain the 
G-type AFM character of interatomic correlations in LaTiOa. 
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Appendix A. Random-Phase Approximation for the Magnetic Transition 
Temperature in Many-Atomic Case 

In this appendix, we present the generalization of the well-known expression for the 
magnetic transition (Curie or Neel) temperature in the random phase approximation 
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(RPA)i 

where J(q) is the Fourier image of Jij and Q is the vector describing the magnetic 
structure in the ground state, to the case of the complex lattices containing n magnetic 
atoms in the primitive cell. 

First, we transform the parameters of the Heisenberg model to the local coordinate 
frame: Jij — Jij = SijJij, where Sij= +1 and —1 corresponds to the FM and AFM 
arrangement of spins in the bond i-j in the magnetic ground state. For a coUinear 
magnetic configuration, this procedure is equivalent to the shift of the origin of the 
Brillouin zone in the right-hand side of ( lA.ll) . Then, we construct the dynamical matrix: 

5(q) = \\JeSu' - -^«'(q)||, 

where 

is the Fourier image of Jij acting between the atomic sublattices £ and M is the 
number of the primitive cells, 

n 

&=1 

and Rj is the radius- vector of the site i. Then, by diagonalizing 5(q), one can find its 
eigenvalues (^^(q) (the "magnon energies") and by repeating the RPA arguments 
derive the following expression for the magnetic transition temperature: 

n^(5 + l) 



2>kBS^ 

Appendix B. Mean-Field Approximation for the Magnetic Transition 
Temperature in the Case of the Noncollinear Spin Arrangement 

The mean- field approximation for the relative magnetization crj(T), which is the 
temperature average of ct,, is formulated in the following way. The mean field (or 
the molecular field), corresponding to the spin Hamiltonian ([2]), is given by 



hi = ^^AijfTj^T). 



Then, the temperature average of &i in the molecular field hi has the following form 
<T) = ^Bs ( ^ (B.i; 



1^,1 \kBT 



® also known as the Tyablikov approximation, renormalized spin-wave theory, or the spherical 
Heisenberg model HOI HZl EH] • 
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where Bs is the Brillouin function for the spin 5*= 1/2 [28j. The equation (IB.ip is solved 
self-consistently. The transition temperature is defined as the minimal temperature for 
which (Ti{T) = 0. 
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